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Abstract 

| In this paper we consider the problem of efficiently constructing fc-vertex fault-tolerant geometric t- 

spanners in R d (for fc > and t > 1). Vertex fault-tolerant spanners were introduced by Levcopoulus et. al 
in 1998. For k = 0, we present an 0(n log n) method using the algebraic computation tree model to find a 
i-spanner with degree bound O(l) and weight 0(uj(M ST)). This resolves an open problem. For k > 1, we 
present an efficient method that, given n points in M d , constructs fc-vertex fault-tolerant i-spanners with the 
maximum degree bound O(k) and weight bound 0(k 2 uj(M ST)) in time 0(n log n). Our method achieves 
the best possible bounds on degree, total edge length, and the time complexity, and solves the open problem 
of efficient construction of (fault-tolerant) i-spanners in R d in time 0(n log n). 
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1 Introduction 



In this work we consider the problem of constructing spanner graphs to approximate the complete Euclidean 
graph. Given an edge weighted graph G = (V, E, W), where w(e) is the weight of an edge e, let do(u,v) 
denote the shortest distance from vertex u to vertex v in graph G. The weight of the graph G, oj(G), is the sum 
of the edge weights of edges in G. A subgraph H = (V,E'), where E' C E, is called a i-spanner of the graph 
' G, if for any pair of vertices u, v G V, dn(u, v ) < t • do(u, v). The minimum t such that H is a t-spanner of G 
■ is called the stretch factor of H with respect to G. An Euclidean graph is a graph where the vertices are points 
in and the weight of every edge (u,v) is the Euclidean distance \\uv\\ between its end-vertices u and v. 
Spanner graphs in R rf have been extremely well studied. We consider spanner graphs with additional properties 
of bounded degree, low weight and fault tolerance. 

In this paper, we study t-spanners and k vertex fault-tolerant t-spanners ((fe, t)-VFTS for short) for a set 
V of n points in M, d . A subgraph H = (V, E) is k vertex fault-tolerant, or fc-VFT for short, if for any pair of 
vertices u and v with uv E, there are k + 1 vertex disjoint paths from u to v in H. Here two paths IIi and II2 
from u to v are said to be vertex disjoint if the only common vertices of IIi and II2 are u and v. A geometric 
graph H = (V, E) is termed (k, t)-VFTS if for any subset F C V of at most k vertices and any two vertices 
w\, W2 G V \F, the graph H(V \ F, E' H ), where E' H = Eh \ {(u, v) \ u G F, or v G F}, contains a path 
U.(wi, W2) from w\ to W2 with length at most 1 1| iwiit»2 1| ■ Given an Euclidean graph that is k vertex fault-tolerant 
(VFT) and a real number t > 1, the aim here is to construct a subgraph H which is a (fc, £)-VFTS subgraph, 
with a bounded vertex degree, and a bounded weight, i.e., uj(H)/uj(MST(G)) is bounded by a specified small 
constant, where MST(G) is the minimum weighted spanning tree of G. 

A greedy algorithm has been used to construct spanners for various graphs lf8l [13^4X71 1271 . For a graph 

G = (V,E) with |V| = n and an arbitrary edge weight, Peleg and Schaffer ll25l showed that any i-spanner 

11 L_ 

needs at least n *+ 2 edges; thus there is edge weighted graph such that any t-spanner of such a graph has 

1 

weight at least 0,(n t + 2 cj(MST)) (the bound is obtained by letting the weight of each edge be 1). Chandra et 
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al. showed that the greedy algorithm constructs a i-spanner of weight at most (3 + ^)n f - 1 - c ■ w(MST) for 
every t > 1 and any e > 0. Regev |[28l proved that the t-spanner constructed by the greedy algorithm has weight 

at most 2e 2 In n ■ -w(MST) when t e [3,2 log n + 1], and has weight at most (1 + 4 lo f f ^ | 2 g° g w ) • a; (MST) 
when t > 2 log n + 1, by studying the girth of the constructed i-spanner. Elkin and Peleg |[T4l recently showed 
that for any constant e, A > there exists a constant (3 = /3(e, A) such that for every n-vertex graph G there is 
an efficiently constructible (1 + e, /3)-spanner of size 0(n 1+x ). 

Constructing i-spanners (Ml H MM MM MM M M MM and (k, t)-VFTS mHQHH for Eu- 
clidean graphs has been extensively studied in the literature. For computing t-spanners of O(l) degree and 
0(u;(EMST)) weight, the current best result Q using algebraic computation tree model is a method with time 
complexity 0(n log 2 nj log log n). An 0(n log n) algorithm which uses an algebraic model together with indi- 
rect addressing has been obtained in IIT7TI . While this model is acceptable in practice, the problem of computing 
low weight spanners in the algebraic decision tree model in time 0(n log n) is still open. We resolve this prob- 
lem and extend the techniques introduced in the first part of the paper to allow us to compute the fc-fault-tolerant 
spanners efficiently. 

In this paper we will also consider constructing (k, i)-VFTS for k > 1 for the complete Euclidean graphs on 
n points V in M. d . The problem of constructing (k, i)-VFTS for Euclidean graphs was first introduced in 11201 . 
Using the well-separated pair decomposition [7j, Callahan and Kosaraju showed that a /c-VFT spanner can be 
constructed (1) in 0(n log n + k 2 n) time with 0(k 2 n) edges, or (2) in 0(nklogn) time with O(knlogn) 
edges, or (3) in time 0(nlogn + c k n) with degree 0(c k ) and total edge length 0(c k • w(EMST)). Here the 
constant c is independent of n and k. Later, Lukovszki l22l presented a method to construct a (k, i)-VFTS with 
the asymptotic optimal number of edges 0(kn) in time 0(n log d_1 n + n/cloglogn). Czumaj and Zhao ifTOl 
showed that there are Euclidean graphs such that any (M)-VFTS has weight at least 0(A; 2 )a;(EMST), where 
EMST is the Euclidean minimum spanning tree connecting V. They then proved that one can construct a (k, t)- 
VFTS using a greedy method H for a set V of n nodes, that has maximum degree 0(k) and total edge length 
0(/c 2 w(EMST)). However it is unknown, given arbitrary k, an Euclidean graph and a pair of vertices u and 
v, whether we can determine in polynomial time if there are k + 1 vertex-disjoint paths connecting them and 
each path has a length at most a given value Notice that this problem is NP-hard when we are given 

a graph G with arbitrary weight function. Czumaj and Zhao further presented a method to construct a (k,t)- 
VFTS for Euclidean graphs in time 0(nk log d n + nk 2 log k) such that it has the maximum node degree O(k) 
and total edge length 0(k 2 log n) • w(EMST) for k > 1. Observe that there is a gap between the lower bound 
0(k 2 ) ■ w(EMST) and the achieved upper bound Oik 2 log n) ■ w(EMST) on the total edge length. 

Our Results: The contributions of this paper are as follows. In the first part of the paper, given a set V of 
points (such input points are called nodes hereafter) in M. d and an arbitrary real number t > 1, we present a 
method that runs in time O(nlogn) using the algebraic computation tree model and constructs a t-spanner 
graph whose total edge length is 0(w(EMST)). The hidden constants depend on d and t, or more precisely, 
the number of cones used in our method, which is 0(( jzj) d )- This solves an open question of finding a method 
with time-complexity 0(?ilogn) in the algebraic computation tree model. The main techniques used in our 
methods are listed below. 

1 . We first apply a special well-separated pair decomposition, called bounded-separated pair decomposition 
(BSPD) which is produced using a split-tree partition of input nodes Q. The split-tree partition uses 
boxes that tightly enclose a set of nodes, i.e., each side of the box contains a point from V. In our 
decomposition, we need to ensure that every pair (X, Y) of separated sets of nodes is contained in two, 
almost equal sized boxes, b(X) and b(Y), respectively, where b(X)(b(Y) respectively) contains only the 

'Private communication with M. Smid. 

2 Edges are processed in increasing order of length and an edge (u, v) is added only if H formed by previously added edges does not 
have k + 1 internally node-disjoint paths connecting u and v, each with length at most t||u — v\\. 



2 



node set X(Y respectively). These boxes are termed floating virtual boxes since they can be positioned 
in a number of ways. An important property of the BSPD that we construct is that for every pair of nodes 
sets (X, Y) in the decomposition, the distance between b(X) and b(Y) is not only not too small (these 
conditions are from WSPD), but also not too large, compared with the sizes of boxes containing them 
respectively. 

2. To facilitate the proof that the structure constructed by our method is a t-spanner we use neighborhood 
cones. At every point x we use a cone partition of the space around x by a set of basis vectors. To guide 
the addition of spanner edges we introduce the notion of General-Cone-Direction for a pair of boxes b 
and b'. This notion ensures that every point x contained inside the box b has all the nodes inside b' within 
a collection of cones C, each cone with apex x. The angular span of the cones C is bounded from above 
by some constant (depending on the spanning ratio t), which ensures the spanner property. 

3. To prove that the structure constructed by our method has low-weight, we introduce the empty-cylinder 
property. A set of edges E is said to have the empty-cylinder property if for every edge e G E, we can 
find an empty cylinder (that does not contain any end-nodes of E inside) that uses a segment of e as its 
axis and has radius and height at least some constants factor of the length of e. We prove that a set of 
edges E with empty-cylinder property and empty-region property has a total weight proportional to the 
minimum spanning tree of the set of end-nodes of edge E. 

In the second part of the paper, given V in M. d , t > 1 and a constant integer k > 1, we present a method that 
runs in time 0(n log n) and constructs a /c-vertex fault-tolerant t-spanner graph with following properties: (1) 
the maximum node degree is 0(k), and (2) the total edge length is 0{k 2 )ui (EMST). This achieves an optimal 
weight bound and degree bound of the spanner graph, which is the first such result known in the literature. The 
second part utilizes techniques introduced in the first part. 

The paper is organized as follows. We present our method of constructing i-spanner in time 0(n log n) in 
Section [2 and prove the properties of the structure and study the time complexity of our method in Section [3] 
In SectionUJ we present and study our method of constructing (k, t)-VFTS. We conclude our paper in section 

El 

2 ^-Spanner in R d Using Compressed Split-tree 

In this section, we present an efficient method in the algebraic computation tree model with time complexity 
O(ralogn) to construct a i-spanner for any given set of nodes V in W 1 for any t > 1. Our method is based 
on a variation of the compressed split-tree partition of V: We partition all pairs of nodes using a special well- 
separated pair decomposition based on a variant of split-tree that uses boxes with bounded aspect ratio. 

2.1 Split Tree Partition of a set of Nodes 

We use d(x,y) to denote the distance between points x and y in R. d in the L p -metric for p > 1. We will focus 
on Euclidean distance here. We define our partition of input nodes V G M. d using a Compressed Split-tree, a 
structure first used in ||3TI . Let Xi be the zth dimension in R. d and x = (x\,X2, • • • , xj) be a point in M d . Then 
an orthogonal box b in R. d is {x = (x\, Xi, • • • , %d) \ Li < x\ < Ri}, where Li < Ri, i = 1, 2, • • • , d, are 
given values defining the bounding planes of the box. Given a box b in W 1 we define the following terminology: 

• \b\ is the number of nodes from V contained in the box b. 

• d(6i, 62) is the Euclidean distance between the boxes b\ and &2> i-e., miiXj-g^^g^ \x — y\\. 

• For a box b, denotes the size of box b, i.e. the length, maxKj<^(i?j — Lj), of the longest side of b. 

• The aspect ratio of a box b is defined as the ratio of the longest side-length over the smallest side-length, 
i.e., max 1 < ij < d (i? i - Li)/(Rj - Lj). 



3 



Given a point set, S, we will refer to the smallest orthogonal box enclosing the point set S as b = \J}(S). 
Such a box b is called enclosing-box of a point set S hereafter. Here the enclosing -box b does not necessarily 
have a good aspect ratio. 

Definition 1 (Tight- Virtual Box) Given a box b = [Q(S) we define a box □(£>) as a tight- virtual box if it has 

the following properties: 

1. □(&) 2 6. it contains b inside 

2. longest side of □(&) w exactly "&(b), the longest side of box b. 

3. □(&) /icw an aspect ratio at most a constant f3 <2. 

Given a set V of n (i-dimensional nodes, let <5(V) be the smallest pairwise distance between all pairs of 
nodes in V. We next define a special split-tree similar to the structure defined in 1711311 . 

Definition 2 (Compressed split-tree) A compressed split-tree, termed CT(V), is a rooted tree of d-dimensional 
boxes defined as follows: 

1. Each vertex u in the compressed split-tree is mapped to a d-dimensional box b = M(u), and associated 
with a tight-virtual box □(&). 

2. The root vertex, termed root, of the tree CT(V) is associated with the enclosing-box M(root) = [FJ(V) 
containing all the nodes in V. Associated with this box is a tight-virtual box n(B (root)) which has a 
bounded aspect ratio f3 < 2 enclosing the box M(root). 

3. Each internal vertex u (associated with a box b = M(u) and the tight-virtual box □(&)) in the tree CT(V) 
has two children vertices, ifb contains at least 2 nodes from V. Consider the two boxes, B[,B' 2 , obtained 
by subdividing b into 2 smaller boxes b'i, 1 < i < 2, cutting b by a hyperplane passing through the 
center ofb and perpendicular to the longest side ofb. Shrink B[,\ < i < 2 to obtain minimum sized 
enclosing-box bi, containing the same set of nodes as B[, i.e., each face ofbi contains a node ofV. Let 

= {&i>&2}- With each box bi in 3(b), we associate a tight-virtual box with an aspect ratio 
at most (3 < 2. Then the children vertices of the vertex u are two boxes b\, 62- Additionally, n(b%) and 
□ (62) are disjoint and are contained inside b. 

4. There is a tree edge from b to every bi G 3(b). Notice that neither b± nor 62 is empty of nodes inside. The 
box bfrom which bi is obtained by this procedure is referred to as the father of the bi, denoted as P(6j). A 
box b that contains only one node is called a leaf box. For simplicity of presentation, we assume that any 
leaf box has a size efor sufficiently small < e <C S(V). 

5. The level of a box b is the number (rounds) of subdivisions used to produce bfrom the root box. The level 
of the box M(root) is then 0. If a box b has level j, then each box in 3(b) has level j + 1. 

In Lemma[n we will show that the tight-virtual boxes n(6j), i = 1,2, can be constructed from □(&) in 0(d) 
time 

The tree CT(V) is called a canonical partition split-tree of V. One difference between our structure and 
the split-tree structure used in is that we associate with each box b in CT(V) a tight- virtual box □(&), while 
in (7J different boxes are used. Another major difference is the floating-virtual-boxes to be introduced later. 

Lemma 1 Given a box b and its associated tight-virtual box n(b), we can find the tight-virtual box n(pi) for 
each children box bi G 3(b) in 0(d) time. 

PROOF. Obviously, for the root box Q](V), we can And a tight-virtual box □(ffl(V)) in 0(d) time. We now 
show that given a vertex b = \J}(S) (enclosing-box of some subset S of nodes in V) in CT(V) and a child 
enclosing-box b\ obtained by subdividing b by a hyperplane h, we can construct a tight- virtual box u(b\) 
from the tight- virtual box □(&) efficiently as follows. Let u(b)h be the box, which contains b\, obtained by 
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partitioning □(&) using the hyperplane h. Assume w.l.o.g that b\ is the one located with the same center as 
u{b)h- Box u{b)h can now be shrunk as b s until one of the sides of the shrunk box b s meets a side of b±. Other 
sides of b s which are larger than can be shmnk to meet b\ if the aspect ratio is not below j3. This gives us 
□ (61). The aspect ratio of □(61) is bounded by the aspect ratio of d(b)h if /3 < 2. ED 

For the purpose of constructing a spanner in M. d , we introduce another box, called floating-virtual box, 
associated with each box in the tree CT(V). For box b, let □(P(6))/ l be one of the two boxes that is produced 
by halving the (longest dimension of) tight- virtual box □ (P(&)) and that contains b inside. Since the tight- virtual 
box □(P(6)) has an aspect ratio bounded by j3 < 2, □(P(6))^ has an aspect ratio bounded by j3 < 2 also. 

Definition 3 (Floating- Virtual Box) Consider a compressed split-tree CT(V) for a set of input nodes V. For 
a box b, a box, denoted as (b), is termed as a floating- virtual box associated with the box b if the following 
properties hold: 

1. it includes the tight-virtual box □(&) of the box b inside, 

2. it is contained inside the parent box P(6) ofb, 

3. it has an aspect ratio at most a constant (3 < 2, and 

4. it is contained inside the box □(P(6))/ l , halved from the tight-virtual box □(P(6)). 

It is worth to emphasize that, for a box b, a floating-virtual box to be used by our method is not unique: it 
also depends on some other boxes to be paired with. It is also easy to show that the floating-virtual boxes of two 
disjoint boxes b\ and 62 will be always disjoint because of the property @]in Definition [3] Table Q] summarizes 
some of the notations used in the paper. See Figure [T]for illustration of some concepts defined in this paper. 



Table 1 : Notations and abbreviations used in this paper. 



WSJ 

B(b) 



0(b) 

d(6i,6 2 ) 
Kh,b 2 ) 



the enclosing-box of a set of nodes S 

the two minimum sized boxes produced by halving b and then shrinking the pro- 
duced boxes to be enclosing-boxes 

the tight-virtual box that contains a box b inside and has an aspect ratio < ft < 2. 
a floating-virtual box containing the tight- virtual box □(&) inside (□(&) and 0(6) 

may be same) and has an aspect ratio < (3. Here 0(6) is not unique in our algo- 
rithm: it depends on the box b' to be paired with for defining edges, 
the size of the box b, i.e., the length of the longest side, 
the Euclidean distance between two boxes b\, 62- 

the edge-distance between two boxes b\, 62- This is equal to d(6' 1 , b' 2 ) where b[ and 

b' 2 are the floating- virtual boxes for bounded-separated boxes b\ and 62. 

the parent vertex of a box b in the tree CT(V) 



Observe that the compressed split-tree proposed here is slightly different from the split-tree defined in 0: 
we define the tight- virtual boxes and also associate a tight- virtual box with some floating- virtual boxes: these 
floating-virtual boxes will be determined by a procedure to be described later. Given that a split-tree can be 
constructed in O(nlogn) steps Q, it is easy to show the following theorem (Theorem©. The proof is similar 
to the proof in |7| and is thus omitted. 

Theorem 2 The compressed split-tree CT(V) can be constructed in time O(dnlogn) for a set of nodes V in 

In the rest of the paper, we will mainly focus on the tight- virtual boxes □(&) for all enclosing -boxes b 
produced in the compressed split-tree. For ease of description, we will refer to b and □(&) by b itself. Thus 
(□(&)) will refer to the same box as 0(6). The difference here is that n(6) has an aspect ratio < f3 while the 
aspect ratio of an enclosing-box b could be arbitrarily large. 
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child box b 1 



child box b2 




boxes 



enclosing box b' 



enclosing box b 

Figure 1: An illustration of several concepts denned in this paper. Here for a box b 2 , depending on the pairing 
box, the floating- virtual box for b 2 could be different. When the box b 2 is paired with the box 62 > the floating- 
virtual box 0(62) is shaded as brown in the figure. When the box b' 2 is paired with the box b[, the floating- virtual 
box \s\(b' 2 ) is shaded as green in the figure. 

2.2 Well-Separated Pair Decomposition and Bounded-Separated Pair Decomposition 

Our method of constructing a spanner efficiently will use some decomposition of all pairs of nodes similar to 
well-separated pair decomposition (WSPD) Q. 

Well-Separated Pair Decomposition (WSPD): Recall that, given two sets of points A,BE M. d , a set 
1Z(A, B) = {{A\ , B\}, {A2, B2}, • • • , {A p , Bp}} is called a well-separated realization of the interaction prod- 
uct A <gi B = {{x, y} I x e A,y E B, and x 7^ y} if 

1. Ai C A, and Bi C B for all i E [l,p]. 

2. AiHBi = for all i & [l,p\. 

3. (Ai (g> n (Aj (g> 5j) = for all 1 < i < j < p. 

4. A® B = Uf = i A* ® Bi. 

5. and i?j is well-separated, i.e., the distance d(.Aj,.Bj) > Qm.ax.(-d(Ai), d(Bi)) for some constant g, 
where of a point set X is the radius of the smallest disk containing X, which is of the same order 
of the size of the smallest box U](X) containing X. 

Here p is called the size of the realization 1Z(A, B). 

Given the split-tree CT(V), we say that a well-separated realization of V ® V is a well-separated pair 
decomposition (WSPD) of V based on CT(V), if for any i, Ai and Bi are the sets of nodes contained in some 
enclosing-boxes, U](iti), U{ G CT(V) and Vi € CT(V), respectively. Here we overuse notations a little 

bit, for a box b £ CT(V r ), we also refer via 6 as the subset of nodes from V contained inside b. In Q, it is 
shown that a well-separated pair decomposition based on a split-tree CT(V) can be constructed in linear time 
0(n) when the split-tree CT(V) is given. 

Theorem 3 [7] Given the split-tree CT{V), we can construct a well-separated pair decomposition based on 
CT(V) in linear-time 0(n) and the realization of this WSPD has size 0(n). 

PROOF. For completeness of presentation, we briefly review the proof here. The algorithm itself is recursive. In 
this proof, when we mention a box b, we always mean the tight- virtual box □ (6) corresponding to the enclosing- 
box used in the split-tree. Each vertex of the compressed split-tree CT(V) has 2 children enclosing-boxes, 
containing point sets denoted by Ai,A^. We construct a well-separated realization for A\ <g> A2. 

We now focus on how to construct a realization for a pair of boxes b and b'. If (b,b') is well-separated, 
i.e., d(b,b') > £max(i9(6), #(£/)), then we are done for this pair. Otherwise, for simplicity, we assume that 
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#(&) > #(&'). In this case, let {q | 1 < i < 2} be the 2 children enclosing-boxes of b in the tree CT(V^). 
We then recursively construct a well-separated realizations of q ® b', for each child box a of 6 and then return 
the union of these realizations as the well-separated realization for the pair (b, b'). Based on the compressed 
split-tree CT(V), we will have a WSPD-computation tree T as follows: 

1. The root of T is (bo, bo), where bo is the box corresponding to the root vertex of CT(V); 

2. For each node (b, b ; ) in T, if b = b' , then it has the following children (bi, bj) where 1 < i < j < 2 and 
b\, 62 are children enclosing-boxes of b in the compressed split-tree CT(V). However, observe that here 
(bi,bi) € T for each children box bi of b. 

(a) If d(b, b') > £max(?9(&), $(b')), then it does not have any children since (b, b') is a pair of well- 
separated boxes. 

(b) If d(b, b') < gmax(i?(6), $(b')) and $(b) > -&(b'), then it has 2 children nodes (a, b') where a is a 
tight-virtual child box of b in the tree CT(V). Recall that here q, 1 < i < 2 is produced by halving 
the longest dimension of b and then shrinking the corresponding boxes into the smallest tight-virtual 
boxes with aspect ratio at most (3. 

(c) If d(b, b') < gm&yL('d(b),'&(b')) and ■#(&) < then it has 2 children nodes (6, q) where q is a 
tight- virtual children box of b' in the tree CT(V). 

For simplicity of analysis, for any vertex (b, U) in the WSPD-computation tree T, we always assume that 
fl(b) > ${b'); otherwise, we reorder them and rename them. A careful analysis in Q show that the WSPD- 
computation tree T has size at most 2n • 2 2 + 4n • (4:Qi) d = 0(n) vertices. The theorem then follows. 

Bounded-Separated Pair Decomposition (BSPD): Observe that a pair of boxes (b, b') in WSPD may have 
a distance arbitrarily larger than max( , i9(6), $(#)), especially, the WSPD produced in (7). To produce spanners 
with low weight, we do not want to include arbitrarily long edges in the spanner, unless it is required. To 
capture such a requirement, we propose a new concept, a bounded-separated pair decomposition (BSPD). 
A well-separated pair decomposition (WSPD) based on a compressed split-tree CT(V) is called a bounded- 
separated pair decomposition, if it has the following additional property: each pair of tight-virtual boxes (b, b') 
in this WSPD satisfies the property of bounded-separation of floating-virtual boxes, i.e., there is a pair of 
floating-virtual boxes 62 = 0(6) and b' 2 = \sl(b') and (62, b' 2 ) = (0(6), nas bounded separation. 

Definition 4 (Bounded Separation) A pair of boxes (62, b' 2 ) has the property of bounded-separation if it satis- 
fies the following properties 

1. Almost Equal-size Property: £i$(6 2 ) < $(^2) < £2^(62)- Here we typically choose constants £2 = 
1/ei = 2. The two boxes 62 an d b' 2 are called almost-equal-sized. 

2. Bounded-Separation Property: £1 max(??(&2), ^(^2)) — ^(62)^2) — 62 max^^), $(b' 2 ))for constants 
1 < Qi < £>2 to be specified later. 

Two boxes b and b' present in the bounded-separated pair decomposition are called a pair of bounded- 
separated boxes in CT(V). The choice of the constants Q2 > Qi > depends on the spanning ratio t > 1 
required. The constants Q2 and q\ will be chosen as specified later to ensure the existence of a pair of bounded- 
separated boxes. Observe the fact that a pair of boxes b and b' is in BSPD does not imply that the distance 
between b and b' is in the range [q\ • max( , !?(6), Q2 • max( , !?(6), i}(b'))]; it is possible that the distance 

d(M') > Q 2 ■ max(^(6),i?(6 / )). However, observe that d(P(6),P(6')) < 02 ■ max(i?(P(6)),i?(P(i/))) for the 
parent boxes P(6) and P(b') of b and b'. 

Given the split-tree CT(V), we then briefly discuss how to connect pairs of nodes to form edges in the 
spanner. Our method is based on the concept of potential-edge-boxes. 

Definition 5 (potential-edge-boxes) Two enclosing-boxes b and b' in the compressed split-tree CT(V) (corre- 
sponding to two vertices in CT(V)) are said to be a pair of potential-edge boxes, denoted as (b, b'), if 
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1. The pair of floating-virtual boxes b 2 = 0(6) and b' 2 = 0(6') have the property of bounded-separation; 
the pair of boxes 62 and b' 2 is called the pair of bounded-separated floating-virtual boxes defining the pair 
of potential-edge boxes b and b'. 

2. None of the pairs (b, P(6')), (P(6), V), (P(6), P(6')), is a pair of potential-edge boxes in CT{V). 

For a pair of potential-edge boxes (b, b'), we use (ED (6), 0(6')) to denote the floating-virtual boxes of b and 
b' respectively that define the pair of potential-edge boxes. From the definition of floating-virtual boxes, we 
have the following lemma. 

Lemma 4 For a pair of potential-edge boxes b and b' and the pair of floating-virtual boxes 62 an d b' 2 defining 
them, the floating-virtual box b 2 = \°}(b) 2 b cannot contain another enclosing-box that is disjoint ofb, i.e. 

&"n& = 0(6) n b" = 0. 

Definition 6 (edge-distance of potential-edge boxes) For a pair of potential-edge boxes b\ and b[, define its 
edge-distance, denoted as £(bi,b'i), as the distance between the pair of bounded-separated floating-virtual 
boxes 62 and b 2 that define the pair of potential-edge boxes b± and b' v i.e., £(bi, b^) = d(b 2 , b' 2 ). 

Given a compressed split-tree CT(V) and associated boxes at each node we then define the edge-neighboring 
boxes of an enclosing-box b as 

Af(b) = {b' I boxes b and b' are a pair of potential-edge boxes in CT(V)} 

Observe that here the distance d(b, b') could be arbitrarily larger than the maximum size of boxes b and b'. 
However, the distance d(P(6), P(6')) < Q 2 max(#(P(&)), tf(P(6'))) if b' £ M{b). Given a fixed size L, and a 
tight-virtual box b, the number of tight-virtual boxes b' such that (1) (\°i(b), \si(b')) has the property of bounded- 
separation, and (2) \s}(b) has size L, is clearly bounded by a constant. However, this does not mean that the 
number of tight- virtual boxes, the cardinality of Af(b), that could pair with b to form a pair of potential-edge 
boxes is bounded by a constant. The reason is that the floating-virtual boxes for a tight-virtual box b depend on 
with which box the box b will be paired (see FigureQ]for illustration). There is an example of nodes' placement 
such that the cardinality |A/"(6)| could be as large as 0(n). The example is as follows: in 2D, we place a node 
at vq = (0, 0) and n nodes V{ at (0, 2 % — e), for 1 < i < n. An additional node v n+ \ is placed at (0, —2"). Here 
e > is a sufficiently small number. Then the smallest box b containing node v n+ \ will have 0(ra) boxes in 
M{b) (the sizes of these boxes are about 2 l , < i < n — 1). 
Thus, our construction method will use another set instead 

M>(b) = {b' I b' e M{b) and #(P(6')) > tf(P(6))} 

Observe that for any pair of potential-edge boxes b and b', we either have b' E Af>(b) or b G M>(b'), or 
both. Consider P(6) and P(6'). The distance between P(6) and P(6') is at most q 2 max(??(62), ^(6' 2 )) since the 
floating-virtual box for b is always contained inside P(6). Thus there are at most ©((^2 • ^2)^) boxes b'. 

Lemma 5 The cardinality ofAf>(b) is bounded by a constant Q((p2 ■ £2)^)- 

We now show that we can construct a linear size BSPD based on CT(V) in linear time. 

Lemma 6 Given any pair of tight-virtual boxes b and b' in a BSPD with $(b') < $(b) < "!?(P(6')), in time 
0(d), we can find a floating-virtual box b" inside F(b')h of almost-equal-size with b. 
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PROOF. Notice that the size of F(b')h is at least $(b)/2 since i?(6) < #(P(b')). Pick any dimension, say d\. 
Projecting F(b')h on this dimension results in a segment, say xy. Let the segment ab be the projection of b' in the 
dimension of d\. Note that [a, b] C [x, y\. We then align a floating-virtual box b" ( s.t. i9(6)/2 < #(&") < #(&)) 
such that its projection on the dimension that contains [c, d] starts at x. If [c, ci] contains [a, b] we are done. 
Otherwise, align b" such that d = b. Since b — x > d — c and d — c > b — a, the alignment of the box 6" 
in dimension di is possible. This can be repeated for all dimensions. It is easy to show that the size of the 
floating-virtual box b" is at most i?(6), and at least i?(6)/2. This finishes the proof. ED 

Theorem 7 Given the compressed split-tree CT(V), we can construct a bounded-separated pair decomposi- 
tion (BSPD) (using constants Q\ and Q2) in linear-time 0(n) and the realization has size 0(n). The constants 
Qi and Q2 are related to q of the WSPD as follows: 



PROOF. We will prove this based on the well-separated pair decomposition computed in the proof of Theorem 
[3] In the proof, we will mainly focus on the tight-virtual boxes, instead of the actual enclosing-boxes. Observe 
that in the WSPD computation tree T (defined in proof of Theorem [3} of a WSPD based on the compressed 
split-tree CT(V), all the leaf vertices will form a WSPD. We first build a WSPD with a constant q, where the 
exact value of q will be determined later. 

Consider a pair of boxes (6, b') at a leaf vertex in the WSPD-computation tree T, i.e. (b, b') is an element of 
the computed WSPD. We show that, by adjusting the sizes of pairs of boxes in the WSPD computed, we can get 
a BSPD. More specifically, for each pair of boxes b and b' in WSPD, we show how to obtain two almost-equal- 
sized floating-virtual boxes b 2 (containing b) and b' 2 (containing b') such that ^1^(62) < d(p2, b' 2 ) < £?2$(^2) 
for some constants Q\ < Q2- 

Assume w.l.o.g., #(£/) < #(&). First of all, because of the properties of the WSPD computation-tree T, 
#(&) < #(P(6')) when (b,b') is a leaf node in the computation tree T. Notice that here, to get the vertex 
(6, b'), we could have split P(ft') first or we could have split F(b) first in the WSPD computation-tree T. Thus, 
by Lemma[6l we can find a floating-virtual box b" = \si(b'), which is almost-equal-sized to b (i.e., i9(6)/2 < 
-d(b") < $(&)), is inside P(6'), and contains b' inside. 

We now show that the distance d(6, b") is at least a constant fraction of $(6). Obviously, 

d(b,b") > d(b,t/)-Vd#(b") > ^max(i?(6),??(6 / )) - Vd#(b) = (q - Vd)ti(b) > Q^(b). 

On the relations of d(6, b") and the size ??(&), there are two complementary cases here: 

Case 1: d(6, b") < Q2"&{b)'- In this case, we have already found a pair of floating-virtual boxes b 2 = b and 
b' 2 = b" for the pair of boxes b and b' such that the distance between the floating-virtual boxes satisfies 
that Qi$(b) < d(b, b") < Q 2 &(b). Thus, we put (b, b') into BSPD. 

Case 2: d(b,b") > Q 2 ${b): Here, from i9(b) < tf(P(&')). we have d(6,P(6')) < £>max(tf(&), tf(P(&'))) = 
Q'd(P(b')) (If this is not true, clearly in the WSPD computation tree T, we will use (b, P(6')) instead of 
(&,&')). Similarly, we have d(P(6), b') < ^max(^(P(6)), 

We now show how to find equal-sized floating-virtual boxes &2 inside the tight-virtual box P(6) and b' 2 
inside the tight-virtual box P(6'). This will identify the potential edge-boxes. Let A be the size of the 
equal-sized boxes 62 and b' 2 . Then we have 




(1) 



d(6 2 ,6 2 ) > d(b, b") - 2VdA 
d(b 2 ,b' 2 )<d(b,b") 



(2) 
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Then the following is clearly a sufficient condition for the existence of such a pair of floating- virtual boxes 
62 and b 2 to define the potential-edge-boxes b and b': 



(d(b 2 ,b' 2 ) > d(b,b") - 2VdA > 6l A 
\d(b 2 ,b' 2 )<d(b,b") < 02 A 

Notice that ■#(&) < d ^' 6 \ Thus, it is equivalent to require A in range 



(3) 



d(b,b") <A< d(b,b") 



92 qi + 2\/d 

Clearly, we have a solution for A when 



Q2 > Qi + 2Vd (4) 

We now show that the boxes 62 and b' 2 will be inside the boxes □ (P(6)/ i and □(P(6')/ l respectively indeed. 
To ensure this, we only need the condition that size A is at most min(-i?(P(&))/2, $(P(&'))/2). 

If tf(P(6)) > t?(P(6')), then we have 

d(b,b") < d(b,¥(b')) + 2Vd$(F(b')) < (q + 2y/d)ti(F(b')) = (2g + 4Vd) min(7?(P(6))/2, tf(P(7/))/2) (5) 
If i?(P(6)) < t?(P(6')), then we have 

d(b,b") < d(P(b),6") + 2\/di?(P(6)) < (g + 2\/d)i?(P(6)) < (2g + 4Vd) min(7?(P(6))/2, tf(P(7/))/2) (6) 
Thus, when 

2^ + Wd < 02 (V) 

we can choose A = < min(7?(P(6))/2,^(P(6 / ))/2) as the final size of 62 and b 2 to ensure that 

the floating-virtual boxes 62 and 6' 2 will be inside the boxes □(P(6))/ l and □(P(6')/ l respectively indeed. 
Observe that we do not put any condition on the locations of 62 and b' 2 here. Thus, as in Lemma [6] we can 
find arbitrary locations of &2 and b 2 in time 0(d). 

Then a sufficient condition for the constants £1 and 02 such that we can compute a BSPD from the WSPD with 
a constant is 

{01 < p — \/d 
r~, (8) 
02 > 20 + 4:Vd 

It is easy to show that we need 02 > 2 01 + 6Vd. This finishes the proof. ED 

Based on the above proof, we can find the floating- virtual boxes for each pair of boxes in the bounded- 
separated pair decomposition in time 0(d). Consequently, we also can find JV> (b) for all boxes b in linear-time. 



2.3 Cones and Cone Partition 

We consider points in the d-dimensional space R d . Let B = {zi, Z2 ■ ■ ■ z g } be a set of g linearly independent 
vectors. The set of vectors B is called a basis of M. d for a cone partition. We define the cone of B, C(B), as 

C(B) = { X iZi I Vt, Xi > 0} 
i<i<g 
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For a point x, we define the cone with apex at x generated by a basis B as C(x, B) = {y | y — x G C(B)}. For 
a given set of points P G R ', the cone region of P is defined as 

C(P,B)= (J C(x,5). 

A vector is said to be in direction of the cone C(f?) (or cone with basis B) if v G C(-u, 5). A vector 
to is said to be in direction of a family of cones with a collection of bases B = {B\,B2, ■ ■ ■ ,B m } if 
v G U^i 1 C(n, 2?j). For two vectors x and y, the angle between them is denoted as <fr(x, y) = arccos( [CTpjbji )> 
where ||x|| is the length of the vector. We define the angular span, based at the origin, of a set of vectors B (and 
its corresponding cone) as 3>(B) = max X]!/ gB{f (i, y)}. Given a set of bases B, its angular span is defined 
similarly. Let J 7 be a finite family of basis of R d . F is called a frame if [j Be jr C(B) = M. d . The angular span 
of a frame F is defined as &(F) = max^gj- §(B). The following lemma has been shown in |[33l . 

Lemma 8 For any < 4> < tt, one can construct a frame F in R d with size 0((^) d ) for a constant c such that 
$(J") < 0. 

The Yao graph, based on the cone partition of the space produced by a frame F, contains all edges uv, 
where v is the closest node to u in some cone C(u, B) for B G F. We use Y(V, F) to denote such graph. The 
following lemma was obvious. 

Lemma 9 The graph Y(V, F) is a t-spanner with 0(n) edges for t = 1 ^2sin(^(J r )/2) w ^ en ^i-^) < 1T /^>- 

It is easy to construct an example of points (for example, n/2 points placed evenly on a side of a unit square 
and another n/2 points placed evenly on the opposite side) such that the weight of Y(V, F) is G(n)w(EMST). 
The maximum node degree in Y(V, F) could also be as large as 0(n), e.g., when n — 1 points evenly placed in 
a circle and 1 node at the circle center. 



General-Cone-Direction Property: We now define the General-Cone-Direction Property for a pair of boxes 
b and b'. Here we require that the pair of boxes b and b' have similar sizes, i.e., e\ ■ #(&) < #(&') < £2 ■ $(b) 
for fixed constants e\, £2 = > 1. In addition, the distance d(6, 6') between the boxes b and 6' is at least 
qi max(i?(6), #(&')) for a constant > 1, where is the size of a box b. We also assumed that the angular 
span of each cone in F is at most a constant a (depending on the spanning ratio t). 

Definition 7 (General-Cone-Direction Property) Given a pair of boxes b and b', a general-cone-direction of 
V with respect to the base box b, denoted as B(b, b'), is a set of cones B C F such that, for every point x G b, 
[J B gB C(x, Bi) properly contains the box b'. Similarly, we can define the general-cone-direction B'(b, b') ofb 
with respect to b'. Also, define 9 as the maximum angular span ofB and B', i.e. 9 = m&x{&(B), <&(B')}. 

A vector uv is said to be in the general-cone-direction of b' with respect to the box b if uv is in the direction 
of B(b,b'). 

Lemma 10 For any pair of bounded-separated boxes b and b' in a Bounded-Separated decomposition, the 
maximum angular span of Bib, b') and B'(b, b') is at most 9 = + 3a. 

PROOF. Figure |2]illustrates our proof that follows. 

A simple computation shows that the general-cone-direction property is satisfied if we select a minimal 
collection B of bases from F such that all vectors in Bi G B are within angle 9\ from vectors in B where 

9 1 <— + a (9) 
Ql 
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Figure 2: An illustration of the proof that we can choose a and q\ such that 9 is small. 



Then 6 < 26>i + a < ^ + 3a. □ 
Clearly, there is only at most a constant number of cones in B and we can find B in constant time. 



2.4 Method for Low Weight t-Spanner Using Split Trees 

Given a set V of n nodes in U. d and any real number t > 1, we now describe our method to construct a structure 
G = (V, E) in time 0(n log n) such that G is (1) t-spanner, (2) each node in V has a degree 0(1), and (3) the 
total edge length co(G) is 0(w(MST)), where MST is the Euclidean minimum spanning tree over V. 

The basic idea of our method is as follows. We first construct the compressed split-tree CT(V) and associate 
geographic information with CT(V). This includes enclosing-boxes, the tight-virtual boxes, and the floating- 
virtual boxes. By Theorem|7]we can construct a bounded-separated pair decomposition (BSPD) in linear time 
using CT{V). For each box b in CT(V), we will select one node as its representative node, denoted as lZ(b), 
used as a gateway node to connect nodes inside this box to nodes at some other boxes in Af(b). To ensure that 
a node is used as a representative node by at most a constant number of boxes (thus ensuring 0(1) degrees for 
these nodes), we apply the following strategy in selecting the representative nodes of boxes. Each leaf vertex 
will have at least one node from V inside. There are at most n — 1 internal vertices in the compressed split-tree, 
and 2n — 1 vertices in CT(V); thus we need at most 2n representative nodes. We will assign 2 credits to 
each node in V. For each leaf box b of CT{V), we choose a node H(b) inside b and charge 1 credit to the 
chosen node. Using a bottom-up approach, for each internal vertex b, we will select a node H(b) from nodes 
contained inside b that has a non-zero credit. Since each internal vertex has at least 2 children vertices, such a 
representative node can always be found. Thus, we have the following lemma. 

Lemma 11 Each node V{ in V is used at most 2 times as a representative node in CT(V). 

Our algorithm works as follows. Given a pair of potential-edge boxes b\ and b[ (defined by a pair of 
bounded-separated floating-virtual boxes &2 and b' 2 ), and their representative nodes u and v respectively, we add 
an edge uv if: 

1. there is no edge xy already added, where x is inside 62 (it is possible that y b 2 ), xy crossing the 
boundary of 62 (thus y is not inside 62), and xy is in the general-cone-direction of the basis £>(&2j b 2 ). Let 
B be the basis such that the representative node v is contained inside the region C(b2, B). Such an edge 
xy is called crossing edge for box 62 in the direction of B, and 

2. there is no crossing edge zw already added, where z is from b' 2 , zw crosses the boundary of b 2 , and zw is 
in the general-cone-direction of the basis B'(b 2 , b^). 

Algorithm Q] presents our method for constructing a t-spanner in M. d with low- weight, and bounded degree 
property. In Algorithm Q] for each enclosing-box b, each basis Bi G F, and each dimension, we store an edge 
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Algorithm 1 Constructing a i-spanner with low-weight 
1: Define a frame F, with a constant c number of bases Bi, Bi, ■ ■ ■ , B c such that the angular span of any base 

is at most a small angle a. The actual value of the angle a will be given later in proofs. 
2: Build the compressed split-tree CT{V) and a BSPD. With each enclosing-box b in the split-tree CT(V), 
we associate a representative node lZ(b). For each enclosing-box b, we also construct M>(b). 
For each enclosing-box b, each basis Bi G F, and each dimension h, we define an array 
CrossingEdge(b, Bi, h). 

3: Sort the edge-distances (see Definition [6]> between all pairs of potential-edge boxes in increasing order. 

(There are a total of 0(n) pairs of potential-edge boxes, thus, the sorting can be done in time 0(n log n).) 

Ties are broken by the actual Euclidean distance between the representative nodes of the boxes. 
4: for(r = lto^2 b M>(b))do 

5: Select the pair of potential-edge boxes b\ and W 1 with the rth smallest edge-distance. Let bi be the 
floating-virtual box containing 61 and b' 2 be the floating-virtual box containing b' x such that bi and b' 2 are 
a pair of bounded-separated floating-virtual boxes. Let u (and v resp.) be the representative node of box 
b[ (and b\ resp.). 

Let B, B' C J 7 be the collection of bases satisfying the General-Cone-Direction Property w.r.t b\ and b\, 
respectively. We then add an edge uv, only if 

1. V-Bj G B, V/i, such that x is inside the box b\ and y is outside of the floating- virtual box 62. there 
is no "crossing" edge xy in CrossingEdge(b\,Bi, h); and 

2. ~iB[ G V/t, such that 2 is inside the box 6' x and io is outside of the floating-virtual box b' 2 , there 
is no "crossing" edge zw in CrossingEdge(b' l , Bi, h). 

After adding edge uv, \/Bi G ,6, V£^ G £>" and V/i < d, we update the array CrossingEdge(b\, Bi, h) 
and the array CrossingEdge^, B\, h) correspondingly (adding an edge uv to the array if the edge uv 
is added to the spanner). 
6: end for 

7: Let G = (V, E) be the graph constructed. 



xy to array CrossingEdge{b,Bi,h) (if there is any) such that (1) node x is inside the enclosing-box b, (2) 
node y is in the cone C{x,Bi), and (3) ?/ is the node that is furthest from the box b in the dimension h if there 
are multiple edges satisfying the first two conditions. This will ensure the following lemma: 

Lemma 12 For every direction specified by the basis Bi, there exists an edge W1W2 with w\ G b and W1W2 
crossing a floating-virtual box 62 (at some ancestor of b) in the direction Bi if and only if there is an edge 
xy G CrossingEdge(b, Bi, h)for a dimension h and xy crosses the virtual box bi (i.e. y bi). 

3 Properties: Low- Weight, Spanner, and Low-Degree 

We next show that the constructed structure G by Algorithm [T] is a t-spanner, has a bounded degree, and has 
low-weight (by choosing the angular span of the frame F , 9, and the parameters Q\, qi, e\, and e 2 in bounded- 
separateness carefully). 

3.1 Degree Property 

Theorem 13 Each node in the constructed graph G by Algorithm^ v G V, has degree < \F\ = 0((^) d ) 
where a is the angular span of the frame F. 
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PROOF. Since each node will serve as a representative for at most two different enclosing-boxes, it suffices to 
show that for each enclosing -box 61, we will add at most a constant number of edges for the representative node 
of this box. We will show that we add at most 1 edge to a node u in any cone direction when u is a representative 
node of a box b\. Assume that we have already added an edge uv in a direction B, where u is the representative 
node in box b\ and v is the representative node of b\ such that b\ and b\ is a pair of potential-edge boxes defined 
by a pair of bounded-separated boxes b 2 (containing 61) and b 2 (containing b^). We show that we cannot add 
another edge uw in the same direction B later. Assume that we did add another edge uw later, because of the 
existence of a pair of bounded-separated virtual boxes ^3 and b' 3 that defines a pair of potential-edge boxes b\ 
and b' (for b' / b^). Then there are only two complementary cases: 

1. 63 contains boxes 61 and f 1 inside. This violates the condition (condition 3) of the potential-edge defini- 
tion: 63 will contain the parent box (which is an enclosing-box) of b\ and b\ inside. Notice that since a 
virtual box 63 contains both b\ and V x inside, our compressed split-tree construction shows that the parent 
box of the boxes 61 and b' x is inside 63 (may be same as 63). 

2. 63 does not contain b\ inside. Then edge uv will be a crossing edge that crosses the boundary of virtual 
box 63. Thus edge uw will not be added. 

This finishes the proof. □ 

Note that we later will show that the angular span a of J 7 depends on the spanning ratio t > 1 that is 
required. 



3.2 The Spanner Property 

We now prove that the final structure G is a t-spanner, where t > 1 is a given constant, if we choose 6, g\ and 
Q2 carefully. 

Theorem 14 The final structure G constructed by Algorithm\J\is a t-spanner for a given constant t > 1 if 'we 
carefully choose a = ^(J 7 ), Q\ and Q2 according to conditions illustrated in f liOl ). 

' qi > tyd 
iVd^ < pi 

< 2^ + (1 + 2^)(1 + ^)/(l - 2sin(0/2) - 2^) < t (10) 
< ^ + 3a 
.02 > 2ei + 6\/d 

PROOF. Note that the last condition is to ensure that we can construct a BSPD. We then prove the theorem by 
induction on the rank of the Euclidean distance between all pairs of nodes u and v from V. First, for the pair of 
nodes u and v with the smallest distance, edge uv clearly will be added to G. Thus, we have a path in G with 
length at most td(u, v) to connect u and v. Assume that the statement is true for all pairs of nodes with the first 
r smallest pairwise distance. Consider a pair of nodes u and v with (r + l)th smallest distance. 

Since we produce a BSPD for the set V of nodes using CT{V), in the box tree CT(V), there will be a 
pair of floating- virtual boxes 62 (containing u) and 6 2 (containing v) that is a bounded-separated pair. Let b\ be 
the largest enclosing-box (from tree CT(V)) that is contained inside b 2 and contains u; and b[ be the largest 
enclosing-box (from tree CT(V)) that is contained inside b' 2 and contains v. Then the pair of boxes 61 and b' x is 
a pair of potential-edge boxes. Depending on whether we have a crossing edge xy when processing the pair of 
potential-edge-boxes b\ and b\, we have the following two complementary cases. 

Case 1: We have an edge xy where x is a representative node of box b\ and y is a representative node of 
box b' x . In this case, we have d(u, v) > d(&i, U-f) > dfo, b' 2 ) > Qi m&x($(b 2 ),$(b' 2 )). By choosing 

Qi>tVd (11) 
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we have d(u,x) < $(61) ■ Vd < d(b2,b' 2 )Vd/ Q\ = £(61, 6' 1 )v / d/£i < d(u,v)\fd/ q\ < d(u,v)/t. Then 
by induction, we have a path connecting x and u with length at most td(u, x): this is true because this path 
can only use edges with length smaller than d(u, v), and d(n, x) < d(u, v ). Similarly, we have a path in G 
connecting v and y with length at most td(v, y). Thus, in the final structure G, we have a path (with subpath 
from u to x, subpath from y to v, and edge xy) connecting u and v with length at most 

td(u, x) + td(v, y) + d(x, y) < (t + l)d(u, x) + (t + l)d(v, y) + d(u, v) < ((2(t + 1) Vd/ft) + l)d(«, t/) 

This is at most id(u, i>) if 

2Vd^j < qi (12) 



P(6i) 



P(6i) 

■ 



Figure 3: An illustration of the proof that the final structure G is a t-spanner. Here for a pair of nodes u and v 
we will have a path with length at most td(u, v), where u and v are representative nodes of the potential-edge 
boxes 61 and 6^. 

Case 2: We do not have an edge xy where x is a representative node of enclosing -box 61 and y is a 
representative node of box b\. In this case, one or both of the following conditions is true: 

1. there is a crossing edge w\W2 such that w\ is inside 61, W2 is outside of floating- virtual box 62 and W2 is 
in the general-cone-direction of C(wi,B), or 

2. there is a crossing edge W1W2 such that W2 is inside 62, w\ is outside of floating- virtual box 61 and w\ is 
in the general-cone-direction of C(w2,B'). 

W.l.o.g., we assume that the first condition is true. See Figure |3]for the illustration of the proof that follows. 

Consider the general-cone-direction £>(6i, 6' x ) of box b[ with respect to the base box 61. The set B of cones 
will be called a meta-cone. The angular span of B is at most a value 8 (from Lemma [TOt. Observe that since the 
meta-cone C(w\,B) will contain the box b\, it will also contain the node v inside. Recall that the edge 101102 
has the same direction as the meta-cone B, the meta-cone C (w\ , B) also contains W2 inside. 

Then for the node W2 and node v, they must be contained in a pair of boxes in BSPD from the definition of 
BSPD. Consider the bounded-separated pair of floating-virtual boxes (say s and s') containing them respectively. 
When the angle < tt/3, we have d(w2, v) < d(w\,v). Together with the fact that W2 is outside of the floating- 
virtual box 62, we can show that the edge-distance (i.e., d(s, s')) of the pair of potential-edge boxes containing 
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W2 and v respectively is less than the edge-distance £(b\, b^) between b\ and b\. In other words, the pair of 
nodes u>2 and v has been processed before the pair of nodes u and v. Thus, we either will have a directed edge 
W3W4 such that W3 and are representative nodes of the boxes s and s' respectively; or we will have an edge 
Z3Z4 such that has the same direction as the meta-cone B(s, s'), i.e., is inside the meta-cone B(s, s'). 
Observe that the distance between the boxes s and s' is smaller than the distance between the boxes b\ and b\ . 

We can repeat the above process and get a sequence of edges w\W2, W3W4, w^we, • • • , W2k-l w 2k' by 
renaming the nodes and the pairs of potential-edge boxes, and the pairs of bounded-separated floating-virtual 
boxes, with the following properties: 

1 . w\ is inside an enclosing -box b\ and w 2 is outside of floating- virtual box b 2 containing b\ (if it is not, we 
can pick the first one W2i~\W2i such that this property is satisfied); node w 2 is inside an enclosing -box, 
called b\, which is inside a floating- virtual box, called b' 2 . The pair of boxes 62 , b' 2 is a pair of bounded- 
separated floating- virtual boxes. Observe that here the boxes b\, 62 and b' 2 may be different from what 
we called at the beginning of our proof. 

2. In general, for i > 1, node w 2 i~i is inside an enclosing-box b 2 i-i which is inside a floating-virtual box 
b2i\ node W2i is inside an enclosing -box b' 2i _ 1 which is inside a floating- virtual box b' 2i . Here, for i > 1, 
the pair of bounded-separated floating- virtual boxes b 2 i, b' 2i contain the pair of potential-edge boxes 62*- 1, 
b ! 2i _i, which is used to define the edge w^i—iwn, i.e., u^i-i (resp. W2i) is a representative node of the 
enclosing-box b 2 i-\ (resp. & 2 i-i)- Notice that here either the box b' 2i _ 1 or the box b 2 i + i could be the 
larger one between them, although both contain node w 2 i + \. We also have that the node w 2 i+i is inside 
the enclosing-box b' 2i _ t for i > 1, while w 2 i is outside of the floating-virtual box b' 2i , for i > 1. 

3. The angle Zvw 2 i-iw 2 i < 9 for a value 9 in Lemma [TOl 
Thus, we have a path 

U W\W2 <^ ^3^4 • • • W2i-\W2i lU2i+1^2i+2 " " " w 2k-\W 2 k V 

to connect the pair of nodes u and v. Here p *~~> (/ denotes a path constructed recursively to connect nodes p 
and q. By induction, we know that the length of path u w\ is at most td(u, w\ 

similarly the length of the path W2i <^ W2i+i is at most m.ax(ty/d'd(b 2i _ 1 ) , t\fd'd{b 2 i + i)) since either (1) w 2 i 
and u>2i+i are inside b' 2i _ l or (2) W2i and 7^22+1 are inside 621+1- 

Notice that max(i?(6 / 2i _ 1 ), $(&2i-i)) < d ( w 2i-i,w 2 ,) from the definition of potential-edge boxes. Addition- 
ally, max(i9(6 / 2 j_ 1 ), i?(&2j-i)) < £2 min(t?(6 2i _ 1 ), #(&2i-i)) since the floating-virtual boxes 6 2 j-i an< ^ &2i-i 
are required to have similar sizes (within a factor e\ = l/e 2 of each other). Then the total length of the path 

W\W 2 W3W4 ■ ■ ■ W2i-2 W2i~\W2i ■ ■ ■ W2k-\W2k IS at most 

()T d(w2i-l,W2i)) ■ (1 H ) 

Thus, the length of the path u w\w 2 W3W4 • • • u>2i-2 W2i-\W2i ■ ■ ■ W2k-iW2k v is at 
most 

(V d^a-i, t^)) • (1 + + ^ • d(«, t/) (13) 

We then bound the length Yli=i d( w 2i-i,w 2 i). From the general-cone-direction property, when 9 < tt/3, it is 
easy to show that 

d(v,w 2 i-i) - d(v,w 2 i) > (1 - 2sm(9/2))d(w2i-i,W2i). 

Since d(v,w 2i+ i)-d(v,w 2i ) < d(w2i,W2i+i) < N/dmax(#(&2i-i),#(&2i+i)) < ^max(d(w 2 i-i,w 2 i),d(w 2i+ i,w 2 i +2 )), 
we have 

d(v,w 2i ) - d(v,w 2i+1 ) > max(d(w 2i - 1 ,w 2 i),d(w 2i+1 ,w 2 i +2 )). 

Qi 
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Then, we have 

k k-1 

d(u, v) + d(wi,u) > d(w 1 ,v) > ^[d(u, w 2i -i) - d(v, w 2i )] + ^[d(u, w 2i+ i) - d(v, w 2i )] 



i=l i=l 

> (l-2sm(9/2)-2^)Y,d(w 2i - 1 ,w 2i ) 
Consequently, the ratio of the length of the path we found over d(u, v) is at most 

—— + (1 + — — 1 + — )/(l - 2sin(0/2 - 2—) < t (14) 

Q\ Ql Ql Ql 

when 9, and £i are chosen carefully (9 is small enough and qi is large enough). 

This finishes the proof of the spanner property. ED 

It is easy to show that we can carefully choose a = ^(J 7 ), Qi and q 2 that satisfy the conditions in (TTTTb . 
Notice that these conditions are weaker than the conditions required to achieve low weight property, illustrated 
in£L7]>. 



3.3 The Weight Property 

We next show that the weight uj(G) of the graph G constructed is 0(ui(MST)). Our proof technique is based 
on the proofs used in lfTTl[T3Tl . Recall that an edge e is added to graph G when we process a pair of potential- 
edge boxes that are defined by a pair of bounded-separated floating-virtual boxes b and b' £ J\f(b). We then say 
that floating-virtual boxes b, b' and the edge e form a dumbbell (as defined in ifTTI ). For a dumbbell formed by 
edge e = (u, v), for both node u and node v, we associate a cylinder with each node, and call it dumbbell head, 
of suitable size. A dumbbell head is a cylinder of radius <5i||in;|| and height ^H^H with < Si <C S 2 <C 1. 
These dumbbell heads are always contained inside the corresponding floating- virtual boxes. Similar to ifTTI . we 
can group edges of G into g = 0(1) groups E\, E 2 , ■ ■ ■ ,E g such that for edges in each group Ei, we have 
Near-Parallel Property: any pair of edges uiu 2 and viv 2 in a group are nearly parallel, i.e., the angle formed 
by vectors u 2 — ui and v 2 — vi is bounded by a constant 9q. This clearly can be achieved using a partition 
based on cones: we first use a constant number of cones to partition the space IR.^ (where the angular span 
of the cone base is at most 9q). Then each cone defines a group of edges: all the (directed) edges uw 
contained in the direction of this cone. 
Length-Grouping Property: In a group, any two edges have lengths that are either nearly equal or differ by 
more than a large constant factor. This can be achieved by first grouping edges into buckets (the ith 
bucket contains edges with lengths in [d' l+1 L, 5 l L] where L is the length of the longest edge and constant 
5 G (0, 1)). Then form a group as the union of every sth bucket (so the edge lengths from different 
buckets differ by at least 5 s factor). 
Empty-Region Property: In a group, any two edges that have nearly equal length of value x are far apart, i.e., 
the distance between end-nodes of these edges are at least eix for some constant e\. Here ei > could be 
any constant (even larger than 1). This clearly can be done by showing that for each edge uv of length x, 
there are at most O(l) edges that are of similar length and are not far apart (that has at least one end-node 
within distance eix of an end-node of uv). Recall that, in our method, for every added edge uv, we will 
only add at most 1 edge for the pair of potential-edge boxes defining uv and the size of the floating-virtual 
boxes b 2 and b' 2 is at least a constant fraction of the edge length d(n, v ). The virtual boxes b 2 , b' 2 used to 
add an edge uv will be used to define the dumbbells of the edge uv. Recall that the virtual boxes will be 
either disjoint or one is completely contained inside the other. This implies that, given any edge uv, there 
is only a constant number of edges xy that are of similar length and are nearby edge uv. 
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Consequently, we have the following lemma: 



Lemma 15 We can group edges ofG into 0(1) groups such that the edges in each group satisfy the preceding 
properties: near-parallel, length-grouping, and empty-region. 

Here the number of groups produced depending on the values q\, Q2, £\, £2, and f3. Recall that the bounded 
aspect ratio is at most (3 < 2 for all the tight-virtual boxes. However, the aforementioned properties do not 
ensure that the total edge weight of edges in a group is 0{1)oj{SMT; l ) where SMTi is the Steiner minimum tree 
connecting the endpoints of edges in E{. We can construct an example of edges H satisfying the aforementioned 
properties such that the total edge weights could be as large as 0(log n)u(SMTi). To prove that the graph 
produced by our method is low-weighted, we need an additional property: 

Empty-Cylinder Property: for eveiy edge uv and its associated dumbbells, there is a cylinder (with the height 
> rj\d{u, v) and radius at least r/2d(u, v) for some positive constants 771 and 772) with axis using some 
segment of the edge uv such that the cylinder is empty of any end-node of edges in the same group. This 
cylinder is called a protection cylinder of the edge uv. 
Observe that the empty-region property does not imply the empty-cylinder property, and neither does the 

empty-cylinder property implies the empty-region property. 



Lemma 16 By carefully choosing q\, and a ( and thus 6 ), according to conditions illustrated in inequality di5D 
every added edge uv by our Algorithm\7]has the empty-cylinder property. 

PROOF. Assume that uv is added due to the pair of potential-edge boxes b\ and b' x , which is defined by a pair of 
bounded-separated (floating- virtual) boxes 62 and b' 2 . Thus 61 and b' x are contained inside 62 and b' 2 respectively. 
Let B be the base such that v G C(&2, B) and let B be the minimal collection of bases such that for any point p 
inside the box 62, b' 2 G C(p, B), i.e., bases that are in the general-cone-direction B. 

Since uv is added, we know that there is no edge xy crossing 62 with x G 62 and xy is in the general-cone- 
direction B. We will show by a simple contradiction that there is a node p, such that the cone C(p, B) is empty 
of nodes w G" 62 with distance d(p, w) < d(p, b'i). If this is not true, consider all the pairs of nodes x and y 
with x G b2, y ^ &2> V £ C(x, B) and d(x, y) < d(x, b^). Let p, q be the pair with the smallest distance among 
all such pairs of nodes x, y. Then edge pq will exist in the graph G, which contradicts the existence of edge uv. 

Since the cone C(p, B) is empty of nodes, then by choosing a large enough q\, we will have a large empty- 
cylinder at the middle of the segment uv. For example, if we let £1 be four times of the value of q\ that satisfies 
condition ©, i.e., 

16\/d 

6 < 3a + — — (15) 
Qi 

then we have an empty-cylinder near the middle of the segment uv with height almost half of the length d(u, v). 
In other words, if condition (TTBl is satisfied, we have ?]i ~ 1/2, and 772 = Qi/{^Q2)- Recall that here q\ and 
Q2 = 0(01) are constants used to define the bounded-separateness of two almost-equal-sized virtual boxes. ED 

Thus, for any edge uv added by our method, we know that there is a cylinder using a segment wz of uv as 
axis with radius at least 772 ) | tit? [ | for a constant 772, wz has length ||w||/2 and in the center of segment uv. 

Definition 8 (Isolation Property) / fi?1/ A set of edges E is said to satisfy the isolation property if 

1. With every edge e = uv G E can be associated a cylinder C(e) whose axis is a segment of uv, and the 
size of the cylinder is not small, i.e., the height is at least rj\d(u, v) and the radius of the basis is at least 
r/2d(n, v)for some positive constants 771 and 772- 



3 Place n nodes evenly on a line and connect every pair of nodes. Then there is a group of edges produced by the preceding 
partitioning will have a total edge weights of 0(log n)uj(SMTi). 
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2. For every edge e, its associated cylinder C(e) is not intersected by any other edge. 

The following theorem was proved in |[T3l by Das et al. . 

Theorem 17 M3\l If a set of edges E in M. d satisfies the isolation property, then u(E) = 0(1)oj(SMT), where 
SMT is the Steiner minimum spanning tree connecting the endpoints of E. 

Based on this theorem, we then show that the graph G produced by our method is also low- weighted. 
Observe that a group of edges from the graph G, partitioned as previously to satisfy the near-parallel, length- 
grouping, and empty-region properties, may not satisfy the isolation property directly. 

Theorem 18 The set of edges Ei that satisfies empty-region property and empty-cylinder property has a total 
weight at most 0{oo{SMT{)) where SMTi is the Steiner minimum spanning tree that spans the vertices in Ei. 

PROOF. We first use the grouping approach to partition the edges into a constant number of groups Ei, 1 < i < 
g, with each group of edges satisfying the near-parallel, length-grouping, and empty-region properties. It now 
suffices to study the weight of a group Ei. We essentially will show that, for each group Ei of edges produced, 
we can remove some edges such that (1) the total length of all removed edges is bounded by a constant factor 
of the total length of the remaining edges, and (2) the set of the remaining edges satisfies the isolation property. 
If these two statements were proven to be true, the theorem then directly follows. Figure |4]illustrates the proof 
that will follow. 
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Figure 4: Illustration of the proof of the low weight property. Here (a) a long edge may intersect the cylinders 
(dark shaded rectangles) of many shorter edges and (b) the length of the similar sized edges that are intersected 
by an edge is bounded. 

Recall that when we added an edge uv to the graph G, edge uv satisfies the Empty-Cylinder Property. It is 
easy to show that the cylinder associated with uv will not intersect any edge xy (with a much shorter length) 
added before uv and xy G Ei. This can be done by shrinking the protection cylinder by at most a small constant 
factor. 

On the other hand, it is possible that edge uv may intersect cylinders of many edges xy G Ei with shorter 
lengths. See Figure @](a) for an illustration of such case. Here the black shaded regions are protection cylinders. 
Given an edge uv, let I(uv) be the set of edges xy G Ei such that uv intersects the protection cylinder of the 
edge xy. We then process edges uv G Ei, starting from the longest edge, to produce E\ as follows: the longest 
edge uv G Ei is added to E\ and update Ei by removing all edges I(uv) from Ei, we repeat this procedure 
till Ei is empty. Clearly, the set of edges E[ satisfies the isolation property and thus has total weight at most 
u(SMT) = 0(u(EMST)). Observe that, the preceding processing of Ei is only for the proof of the low- 
weight property; we will not remove these edges I(uv) from the constructed structure G. Also observe that all 
edges xy G I(uv), where xy has a length at most 5 s d(u, v), s > 1, must be inside a cylinder with axis uv, 
height almost d(uv), and radius at most rj25 s • d(u, v). 

We then show that the length of edges in I(uv) is at most 0(|| For simplicity, we assume that all edges 
in I(uv) are parallel to uv. The rest of the proof will still hold (with different constants) since the edges in Ei 
are almost parallel. Recall that the edges in I(uv) have the length grouping property and the length of eveiy 
edge is at most \\uv\\. Let h{uv) be edges from I{uv) with length in the range of [8 s ' %+1 &(u, v),5 s ' l d{u, v)], 
where i > is any integer. Here s S> 1 and < S <C 1 are positive constants used in deriving the length 
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grouping property. Let be the total length of edges in Ii(uv). We first show that X\ is at most ed(u, v) for 
a small constant < e < 1. Let x\y\, X2IJ2, •• • , x a y a be edges from I\(uv) such that the projection x\ of 
Xi on edge uv is at the righthand side of the projection y' i _ 1 of node yi-\ on edge uv. Then d{xjyj) < 

(1 — r]i)d(u, v) + 25 s d(u, v) < ed(u, v) for a constant e = 1 — 771 + 25 s < 1 when integer s is chosen large 
enough. Here 7/1 is the constant used to define the ratio of the height of a protection cylinder over the length of 
the edge. 

We then show that we will not have edges wz in I\{uv) such that it will have endpoints w such that its 
projection w' on uv is in the interval [x'j,yj] for some 1 < j < a. Figure [4] (b) illustrates the situation for this 
case. This can be proved by choosing a large constant e\ S> 1 in defining the empty-region property: for any 
edge xjyj, there is no edge wz of similar length such that w is within distance e\d(xj,yj) of Xj or yj. 

We can show that the protection cylinders denned by edges in I\(uv) are disjoint, and these protection 
cylinders are also disjoint from the protection cylinder of the edge uv. Obviously, any edge wz from h{uv) 
cannot have node w ov z falling inside the protection cylinders of the edges in Ij(uv) for j < i. Recall that our 
choices of protection cylinders (their sizes) already ensure that any edge wz from h{uv) cannot intersect the 
cylinders of edges from It{uv) with t < i. Thus, the total length of edges in Ii(uv), denoted as Xi, is at most0 

Here X = d(u, v), and 571 Ylt=i Xt is the total height of the protection cylinders denned by edges in Iq{uv) = 
{uv}, I\{uv), l2{uv), ■ ■ ■ , Ii-i(uv). These protection cylinders are empty of nodes, and also empty of edges 
from Ii(uv). Then it is easy to show by induction that, for any i > 1, Ylt=o %t < wX f° r a constant 77 = 
This finishes the proof. ED 

Thus, by choosing the parameters a, qi, Q2, and 6 satisfying the following conditions 

Qi > Ay/d from Theorem [7] 

Q2 > fiVd + 2qi from Theorem [7] 



> 3a + from Lemma M 



(16) 



81 

Qi > t\fd from Theorem [14] 

2\fd^\ < qi from Theorem [H 

+ (1 + ^)(1 + ^)/(l - 2sin(0/2) - 2&) < t from Theorem M 

our structure is a t-spanner, with bounded degree, and is low-weighted. To satisfy the aforementioned condi- 
tions, it suffices to satisfy the following conditions when t < 3 
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> Qyfd + 2q\ from Theorem [7] 

> 3a + from Lemma [16] 

qi > 2y/d^ from Theorem [H 

2t^A + (i + ^)(1 + ^)/(l - 2sin(0/2) - 2^) < t from Theorem M 



(17) 



It is easy to show that we do have solutions for a, q\, 02 and 6. For example, if we let x = Vd/gi, 

j n 1 • -, 25t+l- % /(25t+l) 2 -160t(t-l) . ... , . 

a = x, and 6 = 19x. By choosing 2\x < 1, we get x = — — ^ — - is a valid solution. 



4 With a small additive value whose total length over all i is bounded by X. This is for the case that we may have an edge xy 
from Ii(uv) such that x is outside of the cylinder using uv as axis and radius proportional to S sl d(u, v), and y is inside this cylinder. 
Obviously, the total length of such edges are at most X = d(it, v). 
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Thus, gltji < ^ = < 25t+j clearly is a solution satisfying the last condition here. This solution also 
satisfies the other conditions in Inequalities (fTTT ). Thus the number of cones produced is proportional to (^) d = 
0((jz^) d ) = 0((j^ I ) d ). Then we have our main theorem: 

Theorem 19 Given a set of nodes V in M. d , Algorithm\Tjconstructs a structure that has low-weight ( O (uj(MST) ) ), 
has a constant bounded degree (0((-) d )), and is a t-spanner in time 0((-) d nlogn), where the constant c de- 
pends on the spanning ratio t and is the number of cones needed in our cone decomposition. 

4 (k, £)-VFTS Spanner Using Compressed Boxtree 

In this section, we show how to extend the previous approach to an efficient method for constructing a k > 1 
fault-tolerant t spanner for any given set of nodes V in IR d , t > 1 and k > 1. 

4.1 The Method 

Our approach follows the construction of f-spanners as in the previous sections. We will assume the construction 
of the compressed split-tree CT(V) and a bounded-separated pair decomposition (BSPD) based on CT(V). As 
has already been proven previously (Lemma [TTj, at each box of CT(V), we can choose 1 representative node 
so that each node V{ G V is used at most 2 times as a representative node of some enclosing-boxes. Since we 
want a £>VFTS, we will choose k + 1 representative nodes for an enclosing-box, if it contains at least k + 1 
nodes inside. For easy of presentation, we define various boxes. 

Definition 9 We call an enclosing-box b a k-box if it contains at least k + 1 nodes inside; otherwise it is called 
a non-A;-box. A k-box is called a leaf- A; -box if it is a k-box and none of its children in the compressed split-tree 
is a k-box. 

For each k-box b, we will choose k + 1 nodes contained inside b as the representative nodes 1Z(b) of the 
box b. We will discuss in detail on how to choose 11(b) for a k-box later. For any box b that contains at most k 
nodes inside, all nodes will serve as the representative nodes 71(b) of this box. 

The basic idea of our method is as follows. Consider a pair of potential-edge boxes b\ and b' x , and their 
k + 1 representative nodes (if both boxes are £>boxes). Here b\ and V l are tight- virtual boxes. Consider the pair 
of bounded-separated floating-virtual boxes b 2 and b' 2 , with 62 contained inside the tight- virtual box P(&i) and 
b' 2 contained inside the tight- virtual box P(6' 1 ). We add an edge uv where u is a representative node of b\ and v 
is a representative node of b\, while the following is true: 

1 . there are less than k + 1 disjoint edges of the form xy where x is from b 2 and xy is in the general-cone- 
direction defined by C(x, B) such that B satisfies the General-Cone-Direction property w.r.t. B. Here B 
is the cone basis such that the representative node v is contained inside the region C(&2, B); and 

2. there are less than k + 1 disjoint edges of the form zw where z is from b' 2 and zw is in a general-cone- 
direction B' w.r.t. B' where B' is the basis such that the representative node u is contained inside the 
region C (l/ 2 , B')\ 

Data Structures Used: In our method, for each enclosing -box b, each dimension h, and each cone basis Bi 
of a frame F, we store a set of at most k + 1 disjoint edges xiyi in an array DisjointCrossingEdge, denoted as 
BCE(b,Bi), such that 

1 . Node xi is inside the enclosing-box b, 

2. Node yi is in the cone C(xi,Bi). 

Another set AllCrossingEdge, denoted as ACE(6, Bi), stores all edges xy with x G b, y b and y G 
C(x, Bi). Our method will ensure that the cardinality of ACE(6, Bi) is at most (k + l) 2 for each cone direction 
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B{. Using AllCrossingEdge array ACE(6, Bi) for all cones Bi G B, we can find the maximum number of 
disjoint edges DCE(&i, B) in the general-cone-direction B, which cross the box b\ (one node inside b\ and one 
node outside of 61). This can be done using a maximum matching in the bipartite graph over two sets of nodes: 
one set is all nodes inside b and the other set is all nodes y b with an edge xy for some node x G b. Let 
X(bi) be the set of end-nodes of edges DCE(&i, B) that is inside b\. Let X(bi) be the set of end-nodes of 
edges DCE(6 / 1 , B') that is inside b^. Clearly, for X{b\) and X(bi), each has size at most k + 1. For any node 
u and a direction B, let deg(u, B) be the number of edges incident onto u in the direction of B. These edges 
were added because of processing some pairs of potential-edge boxes with smaller distance. At each step in the 
algorithm, edges are added based on the degree of the representative nodes and edges in the array, DCE. Our 
detailed method for constructing a (k, i)-VFTS is presented in Algorithm [2] 

4.2 Properties: (k, t)-VFTS Spanner, Low- Weight, and Low-Degree 
Lemma 20 The set AllCrossingEdge ACE(6, Bi) has size at most (k + l) 2 when k > 1. 

PROOF. Observe that, for each node u in the box b, we add at most k + 1 edges in any direction B- L . The 
moment the set ACE(6, Bi) has size (k + l) 2 , we must have at least k + 1 disjoint edges crossing the box b 
in the direction Bi. This follows from the observation that we will not add any more edges in our algorithm if 
DCE(6, Bi) is at least k + 1. □ 

We show that the constructed graph is (A;, t)-VFTS. Recall that a structure is called (k, t)-VFTS if for every 
pair of nodes u and v either (1) edge uv is presented or (2) there exist k + 1 node disjoint paths connecting u 
and v such that the length of each path is at most t\\uv \\. 

Theorem 21 The graph G is a (k, t)-VFTS. 

PROOF. We prove that for every pair of nodes u and v, either edge uv G G or there are k + 1 internally vertex 
disjoint paths IIi, il 2 , • • • , H^+i connecting them and each path has length at most t • \\uv\\. In other words, 
G is a (k, t)-VFTS. We will prove this by induction on the distance between b and b' , the bounded-separated 
floating-virtual boxes containing u and v, respectively. We will refer to this distance as the distance between u 
and v. It is easy to show that the pair of nodes v±, V2 with the shortest distance will be in G. 

Assume that for any 1 < j < r, for the pair of nodes whose distance is the j-th smallest, our statement is 
true. We then consider the pair of nodes u and v with the (r + l)-th smallest distance. 

Consider the boxes b and b' such that u G b and v G b' and b and b' are bounded-separated floating-virtual 
boxes and let b' G C(b, B) for some cone basis B G F. Assume for simplicity that the direction defining B is 
horizontal and b is aligned with the coordinate axis. For notational simplicity, let |6| be the number of nodes 
from V that are inside the box b. 

We have three complementary cases depending on the number of points in b and b'\ (1) \b\ = p < k, (2) 
\b'\ = q < k, and (3) p > k + 1 and q> k + 1. 

First consider the case \b\ = p < k. The case \b'\ = q < k follows the same proof. When we processed 
a box b' from J\f>(b), let q = \b'\. If q < k + 1 — g (where g is the number of edges incident on u in the 
direction of B), then edges are added between the node u and all nodes of b'. Thus, we have an edge uv already. 
Otherwise we know that node u is connected to k + 1 nodes, say w\, u>2, • • • , Wk+i with the condition that the 
distance between u and Wi is no more than uv since we processed the boxes b' in M(b) in increasing order of 
distance to b. Additionally, since we only focus on edges in one specific cone direction, we know the length of 
WiV is also less than that of uv. Then the general proof below applies. 

In general, w.l.o.g., assume that neither u nor v are representative nodes in b and b' . In this case, we must 
have p > k + 1 and q > k + 1. Let the k + 1 representative nodes for b and b' be U = {u\ . . . u^+i} and 
U' = {v\, . . . ffc+i} respectively. The following cases arise 
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Algorithm 2 Constructing a (k, i)-VFTS Spanner with Low-weight 
1: Build the compressed split-tree CT(V). For each enclosing k -box b in the tree CT(V), associate it with 
k + 1 representative nodes 71(b). For each box b, we sort the edge-distances £(b, b') for b' G A/>(6). We 
also sort all the distances {£(b, b') \ b' G JV> (6) , for every box 6} in time 0(n log n). 
We then process pairs of potential-edge boxes b and b' in increasing order of their distance £(b, b ! ). 
Initiate the array ACE(6, B t ) and DCE(6, Bi) as empty. 
2: for(r = lto£ 6 .A/">(6))do 

3: Consider an enclosing -box b\ with p representative nodes ui, U2 ■ ■ ■ u p , 1 < p < k+1, and the enclosing- 
box 6' x G AA(6) with g representative nodes u^, u' 2 ■ ■ ■ u' 1 < q < k + 1 such that the distance b[) 
has rank r among all pairs of potential-edge boxes. Let b 2 = 0(61) and b' 2 = ©(fe'i) be a pair of 
bounded-separated boxes. 

Let B G J 7 be the base such that the majority representative nodes of b\ is contained inside C(&2 ; -B)- 
Let B' G J 7 be the base such that the majority representative nodes of b\ is contained inside C(b\, B'). 
Let £>, 0' C J 7 be the collection of bases satisfying the General-Cone-Direction Property w.r.t B and I?' 
respectively. 

4: Let / be the maximum of the cardinality of DCE(&i , B) and DCE(6' 1 , B'). 

5: Let B) be the sorted list of all nodes {u \ u is inside box b, u G" A(6), and deg(ii, i?) < A;}, in 

increasing order of deg(u, B). Thus, we update Yk(b±, B) and Yk(b' l7 B'). 
6: if (p < k) and < k) then 

7: Add all edges u- L Uj for all pairs of U{ and Uj inside b\. Add all edges u^u'j for all pairs of u' { and n^- 
inside 6^. 

8: for each node Ui in b\ with deg(iij, i?) < k do 

9: We add an edge UiUj to a node u'j in 6^ where u'j has the smallest degree deg(ti^, B') among all 

nodes inside b' x . Update the degree for all nodes and arrays ACE and DCE. 
10: end for 

11: for each node u'- in b' x with deg(n^-, B') < k do 

12: We add an edge Uiu'j to a node Uj in b\ where u,i has the smallest degree deg(uj, B) among all 

nodes inside b\ and UiUj was not added before. Update the degree for all nodes and arrays ACE 
and DCE. 

13: end for 

14: end if 

15: if (p > k + 1) and (g > k + 1) then 

16: Add k + 1 - I edges of the form mu'^ where / = max(| DCE(6i, B)\, | DCE(fe / 1 , B')\). Here n i; 

i < k + 1 — /, are the first k + 1 — 1 nodes in Yfc(&i, -B) and u\, i < k + 1 — /, are the first + 1 — I 

nodes in Yfc(&^, B'). 
17: Update the set DCE(6, B^ and ACE(6, Bj) accordingly. 
18: end if 

19: if (p < k) or (q < k), but not both then 

20: Without loss of generality, we assume that p < k and q> k + 1. 

21: Add all edges UiUj for all pairs of Uj and Uj inside b\. Here 61 contains exactly p nodes inside. 

22: for each representative ui of b\ do 

23: Let g' be the cardinality of DCE^ , B'). 

24: If \Y k (b[,B')\ > min(fc + 1 - deg(uj, B), k + 1 - g'), we add min(A; + 1 - deg(ni, B), k + 1 - g') 

edges from Ui to the first k + 1 — deg(u,, i?) nodes in Y^Q/^ B'); otherwise, we add (1^(6^, -B')| 
edges from i±i to nodes in Yk(b' l7 B 1 ). 

25: Update the array Y k (b 1} B) and F fc (bi, B'), and the set DCE(&i, B[) and ACE(6;, 5?) accordingly. 

26: end for 

27: end if 

28: end for 
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29: Let G = (V, E) be the graph constructed. 



Case 1: For all i < k + 1 there are edges U{V{ connecting m G U and V{ G Note that since 6 and b' 
are bounded-separated, d(u, ttj) < d(u, u), for u, S [/ and d(u, < d(it, u) for i>j G £7'. Thus, by 
induction, there exist k + 1 disjoint paths between u and u,j G £7, and there exist k + 1 disjoint paths 
between u and G £7'. Since Vz, UjUj is part of G, by using Mengers theorem it is easy to see that there 
are k + 1 vertex disjoint paths between u and v in U'. It is also easy to show that the length of each of 
such k + 1 paths is at most td(u, v ). 

Case 2: There is an i, such that UjUj does not exist. If I edges of the form uivi are added where Uj and Vi are 
representatives of b and 6' then one of the following subcases arises: 

(a) There are k + 1 — I edges of the form (x, y) where x £ 6 and y G C(x, B), B in the general cone 
direction as B. Let the nodes inside bbe X\,X2 ■ ■ ■ x m , m = k + 1 — I and the edges satisfying the 
preceding condition be a^j/j. Here is a node inside some other box in the direction of B of the 
box b. Note that distance between yi, 1 < i < m and nodes Vj (by measure of the distance £(b, b') 
for two edge-boxes b containing yi and b' containing Vj) is less than the distance between u and 
vj. We can thus apply induction to show that there are k + 1 vertex disjoint paths between u and 

1 < j < m. And by induction there are + 1 node disjoint paths between y« and Vj,\/j,i. 
The result follows from Mengers theorem that there are k + 1 node disjoint paths connecting the 
k + 1 representatives £7 to A; + 1 representatives £7'. 

(b) A similar result is true when there are edges zy, y G b' and z G C(z,B'), where B' in the same 
general direction as B', where b G C(b', B'). 

This finishes our proof that the structure is /c-fault tolerant. Similar to the proof of Theorem [14] we can show 
that for eveiy pair of nodes u and v, each of the k + 1 disjoint paths found in the preceding constructive proof 
has a length at most t\\u — v\\. Thus, the structure we constructed is a (k, t)-VFTS. □ 

Theorem 22 In the graph G = (V, E) constructed by our method, u(E) = 0(k 2 ) ■ w(MST) when k > 1. 

PROOF. Consider the edges added at every node in CT(V). We group edges into 0{k 2 ) groups and will show 
that the total length of edges in each group is at most 0(cj(MST)). 

Note that for a pair of potential-edge boxes b and b', it is possible that 0(k 2 ) edges are added during the 
construction procedure. This happens when there are p < k + 1 (with p = Q(k)) representatives in one box 
during the procedure. As in Section |3]the edges added to connect representative nodes of pairs of potential-edge 
boxes can be partitioned into 0{k 2 ) groups such that the edges in each group satisfy the properties outlined: the 
near-parallel property, length-grouping property, empty-region property, and empty-cylinder property. 

To partition edges into groups, we assume that the representative nodes U inside an enclosing-box b are 
numbered as u\, U2, • • • , u p , where p < k + 1 and the representative nodes U' inside an enclosing-box b' are 
numbered as v±, V2, ■ ■ ■ , v q , where q < k + 1. We group edges sequentially in increasing order of the distance 
between the pair of the potential-edge boxes b and b'. For simplicity, we will only consider one cone base B and 
all edges added in the direction of B. Since there are only a constant number of cone bases, if the edges added 
in the direction of any cone is at most 0(k 2 )u(MST), the total weight of all edges is still 0(k 2 )uj (MST). For 
each edge uv in the cone direction of B, we will put it into one of the k 2 groups: Eij, 1 < i, j < k. 

Notice that, for a pair of potential-edge boxes b and b', we add edges based on rules defined for three 
different cases based on |6| = p, \b'\ = q: Case 1) \b\ = p < k, Case 2) \b'\ = q < k, and Case 3) p > k + 1 
and q > k + 1. 

In case l,\b\ < k, for each m, we add min(A; + 1, q) edges U{Vj for some nodes Vj G U' . Then the edge UiVj 
is added to group Eij for i < k + 1 and j < k + 1. Notice that we also added p(p — l)/2 < k 2 edges inside the 
box b. Observe that each such added edge (inside the box b) has length at most a small constant fraction of the 
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edge UiVj (added to connect representative nodes of b and b'). We will not add these edges to any group E^j. A 
simple charging (charge the total edge length to one such crossing edge U{Vj) method shows that these omitted 
edges have total edge length at most 0(k 2 ) times the total edge length of edges in a group Eij. We will show 
that uj(Eij) is at most w(MST). 

For case 2, for each we will add min(fc + l,p) edges UjVi for some nodes Uj G U and then we similarly 
add the edge UjVi to group Ej^. Thus, all edges added in these cases are put into different groups, i.e., for any 
group Eij, and any pair of potential-edge boxes b and b', we have at most one edge uv with u from b and v 
from b' . 

The third case is that both boxes b and b' have at least k + 1 nodes inside. Let b 2 and b' 2 be the pair of 
bounded-separated boxes that contain the boxes b and b' respectively. In this case, we will add k + 1 — I edges, 
UiVi, where Uj G U and V{ G V. Recall that we added these Z edges because 

1. there are at most I < k + 1 disjoint edges already leaving the floating- virtual box b 2 in the direction of 
the cone B, and 

2. there are at most I < k + 1 disjoint edges already leaving the floating-virtual box b' 2 in the direction of 
the cone B' . 

Notice that those edges were added when processing a pair of potential-edge boxes with shorter distance. Thus 
they have already been put into some groups (at most 21 different groups, since there are at most 21 such edges). 
Then a newly added edge u,iVi is put into a group that is different from those 21 groups. Notice that this is always 
possible since 21 + k + 1 — I < 2(k + 1) < (k + l) 2 when k > 1. Thus, when we put an edge itjUj into some 
group, we know that there is no edge xy in the same group such that x is inside b (resp. b') and edge xy crossing 
the boundaiy of floating-virtual box 62 (resp. b 2 ) in the direction of B (resp. £>')■ Then similar to Theorem [T8l 
we can prove that each group Eij of edges can be further partitioned into a constant number of subgroups such 
that each subgroup of edges satisfies all properties outlined previously, and thus the total weight of all edges in 
each group Eij is at most 0(w(MST)). This finishes the proof of the theorem. ED 

Theorem 23 Each node v G V has a degree 0(k) in the graph G when k > 1. 

PROOF. Notice that for each node u G V, we add edges incident to it only when (1) it is a node contained inside 
some non-fc-box, or (2) when it is a representative node in some k-box. For simplicity, we will only concentrate 
on edges added in the direction of one cone. 

1. We first study how many edges will be added to a node u when u is inside some non-A;-box. When node 
u is inside a non-/c-box, let b u be the largest non-A;-box that contains u inside. Assume that b u contains 
p < k + 1 nodes inside. According to our method, we will add p — 1 edges to other p — 1 nodes inside 
the box b u , and add mm(k + 1 — g u , q) edges for a potential-edge box b' with q nodes inside, where g u is 
the number of edges incident on u, in the direction of B, and with shorter distances. We have thus added 
at most p — 1 + k + 1 = 0(k) edges in the direction of B when we have processed all enclosing -boxes 
b' in N(b), when u is a node inside some non-/c-boxes. 

2. We then study how many edges uv will be added to u when we process a pair of potential-edge boxes 
b and b' such that b is a A: -box, b' is a non-/c-box, such that u is inside the fc-box b and v is inside the 
non-fc-box b' . 

Let b\ = F(b u ) be the smallest enclosing-box that contains u inside and is a k-box. In other words, box b\ 
is a leaf-fc-box. We adopt a charging argument where we assign credits to nodes inside /c-boxes to account 
for edges added to those nodes. With each node inside the k-box is assigned (k + 1) TYPE- 1 -credits for 
each cone direction (another set of TYPE-2-credits is assigned in the next case). 

Now let us see what will happen when we process the enclosing-box b± = b\. The total free TYPE-1- 
credits of this box b\ required to charge of edges in a certain cone direction is \b\\{k + 1) We will prove 
by induction that 
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Lemma 24 Every k-box b will have at least (k + l)(k + 1) — e2free TYPE-1 -credits where e 2 is the 
number of edges that 

(1) are added during the processing of a pair of potential-edge boxes c and d , where c is a descendant 
k-box ofb and d is a non-k-box in M{c), and 

(2) cross the boundary ofb in the given cone direction when we process this box b and its potential-edge 
boxes to add some edge after its children boxes have been processed. 

Proof of Lemma l24l This is clearly true for all leaf-/c -boxes since it has been assigned (k + 1) TYPE-1- 
credits for each node and a given direction, and it has at least k + 1 nodes inside. In the rest of the proofs, 
when we count the number of edges crossing the boundary of b, we will only count the edges that are 
added during the processing of a pair of potential-edge boxes b and b' for some non-A;-box b'. The edges 
added when b' is a k-box will be studied in the subsequent case. 

Observe that when we process b\ and all boxes b[ £ J\f(b\) where b' x is a non-A:-box, the total number of 
edges added to the node u G b\ in any cone B direction is at most k + 1. Consider the case b\ = b\. Let 
b[, b' 2 , b' 3 , • • • , b\, b' i+l , • • • , b' t be all potential-edge boxes that are non-A; -boxes (with respect to b\) in the 
direction of a given cone B, i.e. ,b'- G C (b\ , B) and \ b'-\ < k. We further assume that d{b'- , b) < d(b'j +1 , b) 
for 1 < j < t — 1. Notice that some of these enclosing-boxes may be inside the box P(&i), while some of 
them may be outside of P(&i). Assume that the first i non-virtual potential-edge boxes b' x , b' 2 , b' 3 , ■ ■ • , b' { 
are inside P(&i), while the rest of potential-edge boxes b' i+1 , ■ • ■ , b' t are outside of the box P(6i). There 
are two cases here: x\ = Yl)=i Wj\> k + 1 and x\ = J2)=i Wj\ < k + 1. 

(a) First consider the case xi = Yl)=i Wj\ — k + 1. Assume that / < i is the smallest index such 
that J2j=i Wj\ — k + 1. We consider the number of total free TYPE-1 -credits we will have for 
the enclosing-box P(6i) in this direction. The total TYPE-1 -credits charged to nodes inside b\ 
for adding edges in this direction when processing enclosing-boxes b' 2 , b ; 3 , ■ ■ ■ , b\ is at most 
ik + l)(k + 1) since we add at most k + 1 edges for the "closest" k + 1 nodes inside b'j, j < f 
and we will not add any edges from nodes inside b'j, j > f, to nodes inside b\. Observe that these 
new enclosing-boxes b[, b' 2 , b' 3 , • • • , b' { will also contribute at least (k + l)x± TYPE-1 -credits to 
the enclosing-box P(&i). Thus, the box P(6i) has at least x\ > k + 1 nodes, each with (k + 1) 
free TYPE-1 -credits left. In other words, our claim holds. Recall that a node gets a free TYPE-1- 
credit when it first becomes a node in some k-box and is thus not charged here. Thus, any node v 
inside boxes b' x , b 2 , 6' 3 , • • • , b' { will not be charged any of its TYPE-1 -credits when we add edges vu 
between b'j and b\, j < i. 

(b) We then consider the case that x\ = Y^j=\ Wj\ < k + 1. Let x 2 = X^=i+i Wj\ be the number 
of nodes inside non-/c-boxes outside of b\. In this case, every node inside b\, b 2 , 6' 3 , • • • , b\ will 
connect k + 1 edges to some nodes inside b\. So the total TYPE-1 -credits charged to nodes inside 
b\ is (k + X)x\. The nodes inside U* =1+i 6j will contribute at most min(x2, k + 1 — x\) ■ (k + 1) 
edges since only each of the closest min(x2) k + 1 — xi) nodes will connect up to k + 1 edges to 
nodes inside b\. Thus, the TYPE-1 -credits left by all nodes inside P(fei) is at least \b\\{k + 1) — 
min(A; + l,x\ + x 2 ) ■ (k + 1), while the TYPE-1 -credits contributed by nodes inside boxes b' x , 
b 2 , b' 3 , ■ ■ ■ , b\ is at least x\{k + 1). Thus, the total free TYPE-1 -credits in box P(&i) is at least 

+ 1) — min(A; + 1 — x\,x-i){k + 1). Observe that, when we process the box P(fei), we already 
have e2 = min(A: + 1 — x%, X2){k + 1) edges crossing the box F(b). These edges are added from 
processing pairs of box b\ and non-A;-boxes b'j, j > i. Thus, our claim holds. 
We can show by induction that Lemma l24l holds for all k -boxes. This finishes the proof of lemma. 

For a A;-box b, assume that it already has e 2 edges crossing the boundary of b in the direction of a cone 
B. Then, by our lemma, we know that it has at least (k + l) 2 — e 2 free TYPE-1 -credits. Assume that 
we process pairs of potential-edge boxes b and b'j, for j < t, where b'j is a non-/c-box. For simplicity, we 
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define variable x\ (as number of nodes from bj that are inside P(6)) and X2 (as number of nodes from b'j 
that are outside P(6)) similarly. We only have to add at most mm(/c + 1, x\ + X2){k + 1) — e2 edges to 
nodes inside b. Similarly, we can show that our statement (Lemma |24} holds for the parent /c-box P(6) 
also. 

Observe that in our algorithm, when we decide to add edges uv for a pair of potential-edge boxes, we 
choose a pair of nodes each has the smallest degree in the corresponding direction. Thus, the degree 
difference among all nodes is at most 1. Since we assigned k + 1 TYPE- 1 -credits to every node for 
adding edges in the case uv when u is from a A;-box, and v is from a non-A>box, the maximum number of 
edges added to a node u in this case will be at most k + 1 in any cone direction. 
3. We then study the total number of edges uv that are added to a node u when u is inside some A; -box b\ 
and v is inside some /c-box b\ and the boxes b\ and b' x are a pair of potential-edge boxes. 
Notice that we add at most k + 1 edges when we processing each pair of potential-edge boxes b and b' 
that are A;-boxes. In this case, we will assign 2 TYPE-2-credits for each node inside a £>box. Recall that 
the TYPE-2-credits are different from the TYPE- 1 -credits in the previous case. We charge a node u a 
TYPE-2-credit if an edge uv is added where u and v are from b and b' respectively. Similar to Lemma 
l24l we can prove that 

Lemma 25 Every k-box b will have at least 2(k + 1) — eifree TYPE-2-credits where ei is the number 
of edges crossing the boundary of b in the given cone direction when we start processing this box b and 
its potential-edge boxes to add some edge after its children boxes have been processed. 
We only added at most 1 edge when a node u is served as a representative node of a k-box b\. The 
statement clearly is true for all leaf-fc-boxes. Then consider a non leaf-fc-box b. If we have a fc-box b' 
such that b' is inside P(6), then after processed b and b', F(b) will have at least 2{k + 1) free TYPE-2- 
credits, where b and b' contributed k + 1 TYPE-2-credits each. Similarly we can show that the lemma is 
true when b' is outside of P(fe). 
Thus the theorem follows. □ 

It is also not difficult to show the following theorem. 

Theorem 26 Algorithm\2\can be implemented to run in time Oikcin + C2nlogn), where C2 = &{{j^i) d ) is 
the number of cones partitioned (which in turn is dependent on the spanning ratio t). 

Proof. For each box b and each direction Bi, we store k + 1 disjoint edges in an array DCE(6, Bi) such that 
the end-nodes are furthest from the box in this direction. These edges uv will be sorted based on the distances 
from v to the box b, where u is a node inside b. For a box b, given the array DCE(&i, Bi), and DCE(&2, Bi) 
where b\ and 6 2 are two children boxes of b, we clearly can update the list for box b in 0(k) time from 2 sorted 
lists from children boxes b\ and 62 as follows. We greedily compare the top elements (the link with the furthest 
node) of two children boxes and get the link uv with v being furthest from b in direction B^, the process is 
repeated till we find k + 1 links. Notice that these newly found k + 1 links surely will be disjoint also. Here if a 
node v is connected to multiple nodes u\ and U2 for DisjointCrossingEdge arrays DCE(&x, B) and DCE(&2, B) 
for both children boxes 61 and 62, we will only pick one link UiV for one of the children boxes and discard the 
other. The total time of such processing is 2 • (k + 1). Since there are 0(n) boxes in total, the total time 
complexity of updating the furthest links list DCE can be done in time 0(kn) for a single cone direction. The 
time complexity then follows from the fact that there are C2 cones. ED 

5 Conclusion 

In this paper, we studied the spanner construction for a set of n points in M. d and also fault-tolerant spanners for 
a set of points in W 1 . Our main contribution is an algorithm that runs in time 0(n log n) to construct a (k,t)- 
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VFTS for Euclidean graph with maximum node degree 0(1 + k), and weight at most 0((1 + k) 2 )u(MST) for 
k > 0. All bounds are asymptotically optimum. It remains an interesting future work to extend the method to 
geodesic distance when we are given n nodes on a surface. 
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